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Abstract 



A study of W-pair production in e e annihilations at LEP is presented, based on 877 W W 
(D 1 candidates corresponding to an integrated luminosity of 57 pb~ x at ^=183 GeV. Assuming 



that the angular distributions of the W-pair production and decay, as well as their branching 
fractions, are described by the Standard Model, the W-pair production cross-section is mea- 
sured to be 15.43±0.61(stat.)±0.26(syst.) pb. Assuming lepton universality and combining 
.£h I with our results from lower centre-of-mass energies, the W branching fraction to hadrons is 
determined to be 67.9±1.2(stat.)±0.5(syst.)%. The number of W-pair candidates and the an- 
c3 ' gular distributions for each final state (qq?i^, qqqq, 1vil'V(i) are used to determine the triple 
gauge boson couplings. After combining these values with our results from lower centre-of-mass 
energies we obtain A/c^O.IIIq;^, A^=O.Ollo!ii an d A=— O.lOloi!) where the errors include 
both statistical and systematic uncertainties and each coupling is determined setting the other 
two couplings to the Standard Model value. The fraction of W bosons produced with a longitu- 
dinal polarisation is measured to be 0.242±0.091(stat.)±0.023(syst.). All these measurements 
are consistent with the Standard Model expectations. 
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1 Introduction 



One of the main motivations to double the LEP energy was to study the W boson properties as 
well as the characteristics of the W-pair production process, such as the total cross-section, the 
angular distributions and the helicity structure. The total cross-section for different final states 
can be used to measure the W decay branching fractions in order to test lepton universality 
and to extract information about the Cabibbo-Kobayashi-Maskawa quark mixing matrix. 

In addition to the t-channel //-exchange, W-pair production in e + e~ annihilations involves the 
triple gauge boson vertices WW7 and WWZ which are present in the Standard Model due 
to its non-Abelian nature. Therefore, all the W-pair production properties can be interpreted 
in terms of the Triple Gauge Couplings (TGCs). Any deviation from the Standard Model 
predictions would be evidence for new physics. 



The most general Lorentz invariant Lagrangian |IJ-|J which describes the triple gauge boson 
interaction involving W bosons has fourteen independent terms, seven describing the WW7 
vertex and seven describing the WWZ vertex. This parameter space is very large, and it is 
not currently possible to measure all fourteen couplings independently. Assuming that the 
Lagrangian satisfies electromagnetic gauge invar iance and charge conjugation as well as parity 
invariance, the number of parameters reduces to five, which can be taken as g\, k z , k 7 , A z 
and A 7 |], [|. In the Standard Model (^=k z =k 7 =1 and A z =A 7 =0. Considerations related 
to SU(2)xU(l) gauge invariance, supported by constraints arising from precise measurements 
on the Z° resonance and lower energy data, suggest the following relations between the five 
couplings 

Ak z = — A/s; 7 tan 2 6 W + Agl, (1) 
A z = A 7 , 

where the A indicates a deviation of the respective quantity from its Standard Model value 
and 9 W is the weak mixing angle. These constraints leave only three independent couplings^, 
Ak 7 , Agl and A(=A 7 =A Z ) which are not significantly restricted ]|,[J] by existing LEP Z° data. 

First TGC studies were made at pp colliders, mainly at the Tevatron by CDF ]7] and DO ||, 
using di-boson production. Previous results from LEP used the data recorded at 161 GeV P,pH|] 
and 172 GeV [11, 1^] with an integrated luminosity of ~ 10 pb _1 at each centre-of-mass energy. 



Here we present an analysis of the OPAL data taken in 1997 at a centre-of-mass energy of 
182.7 GeV. The TGC results are combined with those at lower centre-of-mass energies. 

Anomalous TGCs affect both the total W-pair cross-section and the production angular distri- 
bution. Additionally, the relative contribution of each helicity state would be modified, which 
in turn affects the angular distributions of the W decay products. 

In the limit of small W width and no initial state radiation (ISR), the production and decay 
of W bosons is characterised by five angles. These are conventionally Jl|, [| [L3| taken to be: 

1 In our previous publications we analysed the a parameters which are motivated by specific theoretical 
models and are given by, OiB<j>= Ak 7 - Ag^cos 2 8 W , aw<t>= Ag^cos 2 9 W and ctw=^- In this paper, however, we 
prefer to use parameters which are directly related to the W electromagnetic and weak properties 0,0]. 
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the W~ production polar angle?] 9 W ; the polar and azimuthal angles, 9{ and 0*, of the decay 
fermion, W~— >f, in the W~ rest framef]; and the analogous angles for the W + — >f decay, Q\ and 
02- These five angles contain all the information which can be extracted about the helicity 
structure of W + W~ production and decay. The effects of the finite W width and ISR are not 
large at a centre-of-mass energy of around 183 GeV [[L]]. However the experimental accessibility 
of the different angles in W + W~ events observed in the detector, and therefore the sensitivity 
to the TGCs, depends strongly on the final states produced when the W boson decays. 

In this paper we analyse W-pairs in all possible final states, namely, W + W~^ qq£i^, where 
one W boson decays into hadrons and the other decays into leptonsQ, W + W~— > qqqq, where 
both W bosons decay hadronically, and W + W~— > Iv^'Vp, where both W bosons decay into 
leptons. 

Most of the TGC results of this analysis are obtained for each of the three couplings separately, 
setting the other two couplings to zero which is the Standard Model value. These results will be 
presented in tables and log L curves^. We also perform two-dimensional and three dimensional 
fits, where two or all three TGC parameters are allowed to vary in the fits. The results of these 
fits will be presented by contour plots. 

The following section includes a short presentation of the OPAL data and Monte Carlo sam- 
ples, and in section ||| the selection of our W-pair sample is described. The analysis of the 
W-pair cross-section and the W decay branching fractions is presented in section 4, along with 
the interpretation of this analysis in terms of the TGCs. The TGC analysis using the angu- 
lar distributions is done separately for each channel and is described in sections [5], ^| and [7]. 
The combined TGC results are presented in section |[ Section || includes the summary and 
conclusions of this study. 



2 Data and Monte Carlo models 



A detailed description of the OPAL detector has been presented elsewhere [14,15]. The accepted 
integrated luminosity in 1997, evaluated using small angle Bhabha scattering events observed 
in the silicon tungsten forward calorimeter, is 57.21 ±0.15 ±0.20 pb _1 where the first error is 
statistical and the second systematic. The luminosity-weighted mean centre-of-mass energy for 
the data sample is v / i=182.68±0.05 GeV. Part of the TGC analysis involves the 161 and 172 
GeV data samples taken in 1996 at centre-of-mass energies of 161.33±0.05 and 172.12±0.06 
GeV with corresponding luminosities of 9.89±0.06 and 10.36±0.06 pb~ x . 



The semi-analytic program Gentle [TJJ is used to calculate the W + W cross-section for the 

2 The OPAL right-handed coordinate system is defined such that the origin is at the geometric centre of the 
detector, the z-axis is parallel to, and has positive sense, along the e _ beam direction, 9 is the polar angle with 
respect to z and <f> is the azimuthal angle around z. 

3 The axes of the right-handed coordinate system in the W rest-frame are defined such that z is along the 
parent W flight direction and y is in the direction e~ x W where e~ is the electron beam direction and W is 
the parent W flight direction. 

4 Throughout this paper charge conjugate modes are also included and qq implicitly means any pair of 
quark-antiquark of different (or same) flavour which can be produced in W (Z°) decay. 

5 Throughout this paper, logi denotes negative log-likelihood. 



5 



Standard Model case and also for different values of anomalous TGCs. The calculated Standard 
Model cross-section is 15.72 pb at ^=182.68 GeV using the W mass of M w =80.40 GeV/c 2 
measured at the Tevatron[] [[17]]. 

In the analyses described below, a number of Monte Carlo models are used to provide estimates 
of efficiencies and backgrounds as well as the expected W-pair production and decay angular 
distributions for different TGC values. The majority of the Monte Carlo samples were generated 
at y/s = 183 GeV with M\y = 80.33 GeV/c 2 . All Monte Carlo samples mentioned below, unless 



referred to as "generator level", were passed through the full OPAL simulation program JL| 
and then subjected to the same reconstruction procedure as the data. 

A Monte Carlo sample of signal events for the cross-section analysis was produced by the 
Koralw |l!| generator with W + W~ production diagrams (classQ CC03) according to the Stan- 



dard Model. The Koralw generator has the most accurate simulation of ISR and a complete 
treatment of tau polarisation in qqrz/ T events. This sample is used to estimate the W-pair 
selection efficiency Other Monte Carlo samples used for systematic checks were generated with 
the same set of diagrams using the Excalibur |2(J and grc4f [21| programs. To estimate 



the hadronization systematics, Monte Carlo samples were produced by the Pythia P2[ and 



Herwig generators. 



The TGC studies rely mainly on the Excalibur program which has been used to generate 
samples with anomalous TGCs in order to calculate the TGC dependence of the selection 
efficiency and to obtain the expected angular distributions used in the TGC analysis. 

The separation between signal and background processes is complicated by the interference 
between the CC03 set of W-pair production diagrams and other four-fermion graphs, such as 
e + e~ — > WeT7 e , e + e~ — > Z°e + e~ and e + e~ — * Z°Z°. To study the influence of interference effects 
in the four-fermion final states, the grc4f and Excalibur Monte Carlo generators are used. 
In both cases samples were generated using the full set of interfering four-fermion diagrams. 
These four-fermion samples were compared to samples obtained with the same generator using 
only the CC03 diagrams. Four-fermion samples were generated also with anomalous TGCs. 

Other background sources do not interfere with the signal. The main one, Z / 7 — > qq, including 
higher order QCD diagrams, is simulated using Pythia, with Herwig used as an alternative 
to study possible systematic effects. Other background processes involving two fermions in 



the final state are studied using Koralz [24] for e + e —^fi fi , e + e — >t + t and e + e —*vv, 



and Bhwide [25 1 for e + e~— *-e + e~. Backgrounds from two-photon processes are evaluated using 
Pythia, Herwig, Phojet |2l| and the Vermaseren generator It is assumed that the 

present centre-of-mass energy of 183 GeV is below the threshold for Higgs boson production. 



6 The LEP results for the W mass are not used for the TGC measurement, since they have been obtained 
under the assumption that W pairs are produced according to the Standard Model, whereas W production at 
the Tevatron does not involve the triple gauge vertex. 

7 In this paper, the W pair production diagrams, i.e. i-channel v c exchange and s-channel Z°/7 exchange, 
are referred to as "CC03" , following the notation of m . 
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3 Event selection 



W-pair events decay into fully leptonic (W + W~— > luti'Vp), semi-leptonic (W + W~— > qqlz^) 
and fully hadronic (W + W~— > qqqq) final states with expected branching fractions of 10.5%, 
43.9% and 45.6%, respectively. The W + W~ event selection consists of three distinct parts to 
identify each of these topologies. The selections used at 183 GeV are based upon those used 
in the OPAL analysis at \fs ~ 172 GeV [28|. The three selections are exclusive. Only events 
failing the W + W~— > T.i>il'v# selection are considered as possible W + W~— > qq£z^ candidates and 
only events failing both the W + W~ — » Iv^l'Vy and W + W - ^ qq^£ selections are considered 
as possible W + W~— > qqqq candidates. The fully leptonic and semi-leptonic selections are 
separated into the individual lepton types to test charged current lepton universality. 

The detailed selection algorithms are described below. The efficiency values and accepted 
background cross sections are summarised in Tables A] and 



Event 
selection 


+ 

< z/ e e v e 




+ 


Efficiencies [%] for W+W~ -> 

fJ> + VpTV T TV r TV T 


qqez7 e 


qq/^v 


qqrV T 


qqqq 


e^e V e 


67.3 


0.2 


6.5 


0.0 


0.0 


0.6 


0.0 


0.0 


0.0 


0.0 




2.2 


69.9 


5.5 


0.5 


5.8 


1.2 


0.0 


0.0 


0.0 


0.0 


+ 


12.2 


8.7 


58.2 


0.0 


1.6 


8.9 


0.0 


0.0 


0.0 


0.0 




0.0 


1.2 


0.1 


73.7 


5.7 


0.6 


0.0 


0.0 


0.0 


0.0 




0.1 


3.5 


0.6 


10.8 


56.4 


6.8 


0.0 


0.0 


0.0 


0.0 


+ 

T V T T V r 


0.7 


0.9 


4.9 


1.3 


6.2 


44.1 


0.0 


0.0 


0.0 


0.0 


qqeI7 e 


0.0 


0.0 


0.1 


0.0 


0.0 


0.1 


84.0 


0.1 


5.8 


0.1 


qq>^ 


0.0 


0.0 


0.0 


0.0 


0.1 


0.0 


0.1 


86.2 


4.5 


0.2 


qqru T 


0.0 


0.0 


0.0 


0.0 


0.0 


0.1 


4.5 


4.5 


64.6 


0.6 


qqqq 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.5 


84.6 



Table 1: Efficiency matrix for the 183 GeV event selections determined using Koralw (CC03) 
Monte Carlo events. Each entry represents the percentage of generated events for each decay 
channel which are accepted by the different selections. 



The first six background sources listed in Table [| are different four-fermion final states. Their 
expected cross-sections include both contributions from non-CC03 diagrams and the effects 
of interference. The four-fermion background cross-sections for each final state are calculated 
from the difference between the accepted four-fermion cross-section including all diagrams, and 
the accepted CC03 cross-section. For this determination the grc4f Monte Carlo generator 
is used. The associated systematic uncertainty is estimated by comparing the predictions of 
grc4f and Excalibur. At the current level of statistical precision such interference effects are 
small within our experimental angular acceptance. 
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Accepted background cross-sections (fb) 






Event Selection W+W~ -> 




Source 




qqez/ e 


qq/iz/ M 


qqru T 


qqqq 


qqqq 


± 


2 ± 2 


1 ± 2 


6 ± 4 


220 ± 50 


qqeI7 e 


± 


56 ± 28 


1 ± 1 


53 ± 11 


± 


opft 


± 


65 ± 19 


29 ± 4 


74 ± 7 


72 ± 6 


qqz/77 


± 


± 


± 


12 ± 3 


± 


flvv 


40 ± 23 


± 


± 


± 


± 


CM 


9 ± 2 


1 ± 1 


± 


1 ± 1 


± 


Z°/ 7 qq 


± 


61 ± 17 


28 ± 6 


183 ± 22 


1370 ± 150 


Z°/ 7 -> 


6 ± 1 


2 ± 1 


1 ± 1 


6 ± 2 


± 


Two-photon 


30 ± 30 


13 ± 13 


± 


5 ± 5 


6 ± 6 


Combined 


85 ± 38 


200 ± 40 


60 ± 8 


340 ± 27 


1670 ± 160 



Table 2: Background cross-sections for the 183 GeV W + W~ selections. The values for the 
qqqq, qqez/ e and Jv^l'Vp sources are obtained from the difference between the full four-fermion 
and the corresponding CC03 cross-sections. All errors include both statistical and systematic 
contributions. 

3.1 W + W ~— > Ivgl'vgi event selection 

Fully leptonic W + W~ events are identified as an acoplanar pair of charged leptons with missing 
momentum. The event selection is unchanged from that used in reference [ES| (described in 
detail in f29|). In the \fs ~ 183 GeV data sample, 78 events are selected as W + W~— > Iv^Vgi 
candidates. 

The selection efficiency is estimated to be (78.0 ± 2.3)% where the error includes systematic 
uncertainties. The efficiencies for the individual channels are given in Table [l]. The dominant 
systematic uncertainty on the selection efficiency arises from differences observed when com- 
paring Monte Carlo samples which have different implementations of both initial and final state 
radiation effects and the modelling of tau decays. In addition, systematic errors were included 
to account for data/Monte Carlo disagreement (0.8%), the possibility of events being rejected 
due to off-momentum beam particles which can deposit significant energy in the forward de- 
tectors (0.5%), and the knowledge of the trigger efficiency (0.4%). The expected background 
cross-sections from Standard Model processes are given in Table [| 



3.2 W + W — > qq£is£ event selection 

The W + W~— > qq£z^ selection consists of three separate selections, one for each lepton flavour. 
The W + W - ^ qqT^r selection is applied only to events which fail both the W + W~— > qqez/ e 
and W + W~— > qq//F^ selections. 

The W + W~ — > qq£z^ event selection for the 183 GeV data is a modified version of the 172 GeV 
selection described in detail in |p8|j . A looser set of preselection cuts is used since the lep- 
ton energy spectrum is broader at 183 GeV due to the increased boost. The set of vari- 
ables used in the likelihood selections has been also modified. In the W + W - ^ qqrz7 T sam- 
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pie there is a significant background from hadronic decays of single-W four-fermion events 
(e + e~ — > Wez7 e ) and an additional likelihood selection is used to reduce this. This selection is 
only applied to W + W~— > qqrV T events where the tau is identified as decaying in the single- 
prong hadronic channel. Finally, in order to reduce the e + e~ — > Z°Z° background, events pass- 
ing the W + W~— > qqeI7 e likelihood selection are rejected if there is evidence for a second ener- 
getic electron. A similar procedure is applied to the W + W - ^ qq/ilV selection. 

In total, 361 events are identified as W + W~— > qq£z/£ candidates, of which 140 are selected 
as W + W~— > qqe!7 e , 120 as W + W~— > qq^V^ and 101 as W + W~ — > qqrV T . The efficiencies of 
the W + W~— > qqivg selection for the individual channels are given in Table p]. The efficiencies 
include corrections (~ 0.5%) which account for observed differences between the data and the 



Monte Carlo simulation as described in pq| . The efficiencies are also corrected by a factor 
of 0.991 ± 0.003 for the W + W~— > qqeU e and W + W~— > qq/il 7 ^ selections and by a factor of 
0.985 ± 0.005 for the W + W~— > qqrV T selection to account for events being rejected due to the 
presence of off-momentum beam particles which can deposit significant energy in the forward 
detectors. A systematic error is assigned for possible tracking losses which are not modelled 
by the Monte Carlo simulation of the OPAL detector. This uncertainty is evaluated using 
selections designed to identify W + W~— > qqez7 e and W + W~^ qqT-t^V events where the track 
was either not reconstructed or was poorly measured. Possible biases due to fragmentation 
uncertainties are studied by comparing fully simulated Monte Carlo W + W~^ qqlz^ samples 
with fragmentation performed using Pythia and Herwig. Other systematics are estimated 
using samples generated with different Monte Carlo models and different input W boson masses 
and beam energies. Table |3| lists the sources of the systematic uncertainties on the selection 
efficiencies. 







Signal efficiency error 


(%) 






Event selection W + W 






Source of uncertainty 


qqe^e 




qqru T 


qqqq 


a) 


Statistical 


0.3 


0.3 


0.4 


0.2 


b) 


Comparison of MC models 


0.6 


0.4 


0.5 


0.3 


c) 


Data/Monte Carlo 


0.6 


0.5 


0.6 


0.6 


d) 


Tracking losses 


0.4 


0.4 


0.4 




e) 


Detector occupancy 


0.3 


0.3 


0.3 




f) 


Fragmentation 


0.2 


0.1 


0.6 


0.8 


g) 


M w dependence (±0.09 GeV/c 2 ) 


0.1 


0.2 


0.1 


0.1 


h) 


Beam energy dependence 


0.1 


0.1 


0.1 


0.1 


Combined 


1.1 


0.9 


1.2 


1.0 



Table 3: Sources of uncertainty on the W + W — > qqiug and W + W — > qqqq selection efficien- 
cies. 



Table |2| shows the corrected background cross-sections and their total uncertainties for the 
three qq£z/^ selections. The systematic errors on the expected background cross-sections are 
dominated by differences between data and Monte Carlo for the two-fermion backgrounds and 
by differences between generators in the case of the four-fermion backgrounds. The systematic 
errors on the four-fermion backgrounds were estimated by comparing the expectations of grc4f , 
Excalibur and Pythia. 



9 



An initial estimate of the Z°/7 — > qq background is obtained from Monte Carlo. This estimate 
is then corrected using data to account for possible uncertainties in the Monte Carlo modelling, 
e.g. the lepton fake rate in Z°/7 — > qq events. Several methods are used to estimate the 
correction factors to the Monte Carlo estimates of the Z°/7 — > qq background. The central 
value is obtained using 'fake' Z°/7 — > qq events formed by boosting Z° — > qq events from at 
y/s = 91.2 GeV along the z-axis according to the invariant mass distribution of the hadronic 
system in Z°/7 — > qq events at ^=183 GeV f3*0|| . This procedure is applied to both data 
and Monte Carlo. The correction factor is determined from the relative fractions (data/Monte 
Carlo) of 'fake' events which pass the W + W~— > qq£^ selections. As a result, the Monte Carlo 
estimates of the Z°/7 — > qq background are scaled by a factor of 1.3±0.3 for the W + W~— > qqez/ e 
selection, 1.0 ± 0.1 for the W + W~— > qq/i^ selection and by 1.1 ± 0.1 for the W + W~— > qqrz/ T 
selection. The errors reflect the spread of values obtained using alternative methods, e.g. 
using events passing the preselection and then fitting Monte Carlo signal (W + W~— > qq^) 
and background (Z°/7 — > qq) components to the kinematic distributions of the data. 

In the W + W~^ qq^T selection there is a non-negligible background from hadronic decays of 
single-W events (e + e~ — > Wez7 e ) where the electron goes undetected down the beam pipe. This 
background is determined by the number of fake tau candidates formed out of the fragmentation 
products in the hadronic decay of the W. A correction factor to the Monte Carlo background 
is obtained using fake single-W events formed by removing the lepton candidate from selected 
W + W~^ qqe77 e and W + W~^ qquX/ M events (both in data and Monte Carlo). The ratio of 
W + W~^ qqT77 T selection efficiencies (data/Monte Carlo) for these fake events, 0.8 ± 0.1, is 
used to scale the Monte Carlo estimate of the e + e~ — > Weu e background. 



3.3 W + W — ► qqqq event selection 

The selection of fully hadronic W + W - ^ qqqq events is performed in two stages using a cut- 
based preselection followed by a likelihood selection procedure. The general features of this 
selection are similar to those used previously at 172 GeV f28| , although it has been re-optimised 
to improve the rejection of the dominant background arising from hadronic Z° — > qq decays. 

All events which are classified as hadronic events []3T| and which have not been selected by 
either the W + W~— > Ivfl'Vpj or W + W~— > qq£z/£ selections are considered as candidates for the 
W + W _ — > qqqq selection. Tracks and calorimeter clusters are combined into four jets using 
the Durham k± algorithm 3^] and the total momentum and energy of each of the jets are 



corrected for double-counting of energy |3"3]| . To remove events which are clearly inconsistent 
with a hadronic W + W~^ qqqq decay, candidate events are required to satisfy the following 
preselection criteria: 



The fitted invariant mass of the event scaled by the collision energy, Js'/s, must be 
greater than 0.75. 

The visible energy of the event must be greater than 0.7-^/i. 

The energy of the most energetic electromagnetic cluster must be less than 0.3^/i. 
Each jet is required to contain at least one charged track. 
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The quantity^ j ang must be greater than 0.05. 

The logarithm of the QCD matrix element W420 0) calculated using the jet momenta 
as estimates of the parton momenta, is required to be less than zero. W420 is an event 
weight formed from the tree level 0{a 2 s ) matrix element [|35| for the QCD process (e + e~ — > 
qqqq, qqgg). It is assigned the largest value of any permutation of associated jets to 
partons. 



This preselection rejects 97.7% of the Z° — > qq events which comprise the dominant background 
in the W + W _ — > qqqq channel. The preselection efficiency for the hadronic W + W~— > qqqq 
decays is estimated to be 88.8%. In total, 524 candidates pass the preselection. 

Events satisfying the preselection cuts are classified as signal or background based upon a four 
variable likelihood selection. The following likelihood variables are used since they provide 
a good separation between the hadronic W + W~— > qqqq signal and the dominant Z° — > qq 
background process for a minimum number of variables used: 



log (2/45), the logarithm of the value of the Durham jet resolution parameter at which an 
event is reclassified from four jets to five jets. 



log(W42o), the logarithm of the QCD matrix element |34|. 
I cos#n-r|) the cosine of the modified Nachtmann-Reiter angl 
the sphericity of the event. 



36 



Rather than using the product of the individual probability density functions to construct a 
classic likelihood discriminator, a coordinate transformation technique has been developed to 
reduce correlations between the four input variables [p7|l . 

An event is then selected as a hadronic W + W~^ qqqq decay if the likelihood discriminant 
variable L, shown for data and Monte Carlo in Figure [l], is greater than 0.36. This cut value 
was chosen to maximise the product of signal purity and efficiency. 

In the OPAL 172 GeV cross-section measurement [p8[| , a 5% improvement in the statistical 
error was achieved by counting the event weights of all preselected events rather than simply 
specifying a cut on the discriminant variable. In this analysis, the use of event weights improves 
the statistical error on the cross-section measurement by less than 2% and results in an increased 
systematic uncertainty. For this reason, the more straightforward counting method is used. 

The overall efficiency of the W + W~^ qqqq event selection is estimated from the Koralw 
Monte Carlo simulation to be (84.6±1.0)% where the error is an estimate of all known systematic 

8 With the jets ordered by energy, the quantity j ang is defined as -^=(1 — C12C13C23) where = cos 0jj, 9ij 
being the angle between jets i and j, and E4 is the energy of the fourth (lowest energy) jet. This quantity is 
strongly peaked towards zero for the dominant Z°/7 — > qq background process. 

9 With the jets ordered by energy, the quantity cos6>n_F(. is defined as ^i~_fi^l?^~?^ . This variable, which is 
sensitive to correlations between the underlying parton momenta, tends to be flat for the dominant Z°/7 — > qq 
background and somewhat peaked for the W + W~^ qqqq signal events. 
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uncertainties; the individual components are listed in Table |[ The total expected background 
cross-section from non-CC03 diagrams is estimated to be (1.67 ± 0.16) pb where the error 
represents the systematic uncertainty. The processes contributing to the background cross- 
section are listed in Table |2|. In total, 438 candidates pass the selection. 

The main source of uncertainty for both the signal efficiency and the background cross-section 
is related to the modelling of the fragmentation process. This uncertainty is estimated by 
comparing the selection efficiency for both signal and background events using an alternative 
QCD Monte Carlo model (Herwig). In addition, the parameters cr q , b, Aqcd, and Q of the 



Jetset fragmentation model are varied by one standard deviation about their tuned values |38 
The Monte Carlo modelling of the data is further studied by comparing the distributions of 
the four likelihood variables seen in the data with various Monte Carlo samples. The signal 
efficiency evaluated using Koralw is compared to alternate generators (Excalibur, Pythia 
and grc4f ) to test the Monte Carlo description of the underlying hard process. Uncertainties 
related to the beam energy and W mass are evaluated with Pythia samples generated over a 
range of values. In each case, the observed differences are taken as an estimate of the systematic 
uncertainty. 



4 W + W cross-section and W decay branching fractions 
4.1 Cross-section and branching fraction results 

The observed numbers of selected W + W~ events are used to measure the W + W~ production 
cross-section and the W decay branching fractions to leptons and hadrons. The measured 
cross-section corresponds to that for W pair production via the CC03 diagrams. 



Selected as 


Efficiency [%] 


Purity [% 




Expected 


Observed 


w+w--> ivttpff 


78.0 ±2.3 


93.8 




78.9 ±2.3 


78 


W+W"^ qqtz^ 


84.7 ± 1.3 


90.2 


371.2 ±8.2 


361 


W + W~^ qqqq 


84.6 ± 1.0 


78.2 




442.6 ± 12.4 


438 


Combined 


84.4 ±0.8 


85.2 


892.7 ± 19.2 


877 



Table 4: Observed numbers of candidate events in each W + W~ decay topology for an integrated 
luminosity of (57.21 ± 0.25) pb~ x at (182.68 ± 0.05) GeV, together with expected numbers of 
events, assuming Mw = (80.40 ± 0.09) GeV/c 2 and Standard Model branching fractions. The 
expected numbers of events include systematic uncertainties from the efficiency, luminosity, 
beam energy, W + W~ cross-section (2%) and Mw- In this table cross contamination between 
W + W~ events from different topologies has been included as background (i.e. excluded from 
the efficiency numbers). The errors on the combined numbers account for correlations. 

Table f| summarises the event selections in the three W + W~ decay topologies. The efficiencies 
refer to CC03 W + W~ events. The expected number of events is calculated using the Gentle 
cross-section of 15.72 pb assuming a W mass of 80.40 GeV/c 2 |17| and a centre-of-mass energy 
of 182.68 GeV. The systematic uncertainties on the expected numbers of events include contri- 
butions from the current errors of ±0.09 GeV/c 2 on Mw and ±0.05 GeV on the centre-of-mass 
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energy (both below 0.1%) and a 2% theoretical uncertainty on the cross-section calculation. 
The data are consistent with the Monte Carlo expectation. 

The selected numbers of events can be used to determine the W + W~ CC03 production cross- 
sections separately into each channel: 

a(W + W~^lvd'V e ) = (1.64 ± 0.20 ± 0.07) pb, 
a(W+W~^ qqlu e ) = (6.68 ± 0.39 ± 0.12) pb, 
o-(W+W--* qqqq) = (7.07 ± 0.43 ± 0.21) pb. 



As in |28|| , the W + W~ cross-section and branching fractions are measured using the observed 
events from the ten separate channels summarised in Table |^. Three different maximum like- 
lihood fits are performed. In the first case, <7ww(183 GeV), Br(W — > ez/ e ), Br(W — > jiV^) and 
Br(W — > tV t ) are extracted under the assumption that 

Br(W e77 c ) + Br(W -> fiV^) + Br(W -> rv T ) + Br(W -> qq) = 1. 

In the second fit, charged current lepton universality is imposed and in the third fit, the W + W~ 
cross-section is determined assuming Standard Model branching fractions. The results are 
summarised in Table |6]. 



Selected as 


Expected signal 


Expected back. 


Total 


Observed 




e + ^ e e u e 


8.6 ±0.4 


0.2 ±0.4 


8.8 ±0.6 


12 




e+z/e/i - ^ 


17.5 ±0.6 


0.2 ±0.1 


17.8 ±0.6 


11 


w + w~-> 


e + v c T~V T 


16.7 ±0.6 


1.0 ± 1.0 


17.6 ± 1.6 


20 


w+w--> 




9.3 ±0.3 


1.4 ± 1.0 


10.7 ± 1.0 
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14.6 ±0.6 


1.1 ± 1.0 


15.7 ± 1.2 


15 


w+w-^ 


T + V T T~V T 


7.4 ± 0.4 


0.8 ±0.6 


8.2 ±0.7 


7 


w+w-^ 


qqez7 c 


118.4 ±2.8 


11.3 ±2.3 


129.8 ±3.6 


140 


w+w~-> 


qq/il7 M 


119.8 ±2.7 


3.4 ±0.5 


123.2 ±2.8 


120 


w+w-^ 


qqrz7 T 


99.1 ±3.0 


19.2 ± 1.6 


118.2 ±3.3 


101 




qqqq 


347.1 ±8.3 


95.5 ±9.2 


442.6 ± 12.4 


438 


Combined 


758.4 ± 16.5 


134.3 ±9.6 


892.7 ± 19.2 


877 



Table 5: Observed numbers of candidate events in each W + W~ decay channel for an integrated 
luminosity of (57.21 ± 0.25) pb -1 at (182.68 ± 0.05) GeV together with expected numbers of 
signal and background events, assuming M w =80.40 ± 0.09 GeV/c 2 . The predicted numbers 
of signal events include systematic uncertainties from the efficiency, luminosity, beam energy, 
W + W~ cross-section and M w , while the background estimates include selection and luminosity 
uncertainties. The errors on the combined numbers account for correlations. 



The hadronic branching fraction can be interpreted as a measurement of the sum of the squares 
of the six elements of the CKM mixing matrix, |V;j|, which do not involve the top quark p]: 

Br(qq) 



Br(qq)) 
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Fitted parameter 
183 GeV data 




Fit assumptions : 
Lepton universality 


SM branching fractions 


Br(W -> eV e ) 
Br(W -> 
Br(W -> r77 T ) 
or^W — > tz/gj 
Br(W -> qq) 


0.121 ± .010 ± .003 

A -1 A ^Tf 1 AAA 1 a, a o 

0.107 ± .009 ± .003 
0.094 ± .011 ± .003 

0.678 ± .013 ± .005 


0.108 ± .004 ± .002 
0.676 ±.013 ±.005 




cxww(183 GeV) [pb] 


15.33 ±0.61 ±0.27 


15.43 ±0.61 ±0.26 


15.43 ±0.61 ±0.26 


Fitted parameter 
161-183 GeV data 




Fit assumptions : 
Lepton universality 


Br(W -> ez7 e ) 
Br(W -> ^) 
Br(W -> rz7 T ) 
Br(W -> ?z/*) 
Br(W -> qq) 


0.117 ± .009 ± .002 
0.102 ± .008 ± .002 
0.101 ± .010 ± .003 

0.680 ± .012 ± .005 


0.107 ± .004 ± .002 
0.679 ±.012 ±.005 



Table 6: Summary of cross-section and branching fraction results from the 183 GeV data and 
the branching fraction results from the combination of the 161 GeV, 172 GeV and 183 GeV 
data. The results from three different fits described in the text are shown. The correlations 
between the branching fraction measurements from the fits without the assumption of lepton 
universality are less than 27%. 



where a s (Mw) is taken to be 0.120 ± 0.005. The branching fraction Br(W — > qq) from the 161 
- 183 GeV data obtained from the fit assuming lepton universality yields, 

£ |Vij| 2 = 2.04 ±0.11 ±0.05. 

i=u,c; j'=d,s,b 

This is consistent with a value of 2 which is expected from unitarity. Using the experimental 

^ ^ 2 2 2 2 2 

knowledge [B9] of the sum, |V ud r±|V us r±|V ub r±|V cd | z ±|V cb | z = 1.05±0.01, the above result 
can be interpreted as a measure of |V CS |, which is the least well determined of these elements: 

IVJ = 0.99 ±0.06 ±0.02. 



The measured W + W~ production cross-section at y/s = 182.68 GeV is shown in Figure ^ 
together with the recent OPAL measurements of oww at \/s = 161.3 GeV |4(| and at \/s = 
172.1 GeV [28|. Figure ^ also shows the Gentle prediction which is in excellent agreement 
with the data. On the other hand, the cross-section calculated without the contribution of 
the WWZ vertex, corresponding to anomalous couplings Agl=-1 and Ak z =-1 (dashed line in 
Figure Q) fails to describe the data. 



4.2 TGC analysis using the cross-section 



A quantitative study of TGCs from the W-pair event yield is performed by comparing the num- 
bers of observed events in each of the three event selection channels with the expected number 
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which is parametrised as a second-order polynomial in the TGCs. This parametrisation is based 
on the linear dependence of the triple gauge vertex Lagrangian on the TGCs, corresponding 
to a second-order polynomial dependence of the cross- sect ion. The polynomial coefficients are 
calculated from the expected cross-section in the presence of anomalous couplings determined 
with Gentle and the slight dependence of our selection efficiency on the TGCs is obtained 
from the Excalibur Monte Carlo samples. The Standard Model values for the W branching 
fractions are used. The background, which originates predominantly from Z°/7 — > qq events, 
is assumed to be independent of the TGCs. 

The probability to observe the measured number of candidates, given the expected value, is 
calculated using a Poisson distribution. The product of the three probability distributions 
corresponding to the three event selection channels is taken as the cross-section likelihood 
function. 

The following sources of systematic uncertainty on the expected number of events are consid- 
ered. 



• The theoretical uncertainty in the expected cross-section which is obtained by comparing 
the cross-sections obtained from the Gentle and Excalibur programs. This uncertainty 
depends on the anomalous couplings and has a typical size of 2%. 

• The small effect of the uncertainties in the W mass from the Tevatron measurement, 
(80.40 ± 0.09) GeV/c 2 0, and the LEP centre-of-mass energy, (182.68 ± 0.05) GeV, on 
the total cross-section (less than 0.1% each). 

• The uncertainties in the selection efficiencies and accepted background cross-sections, as 
listed in Tables and £|. 

• The uncertainty in the luminosity, of 0.5%. 



The systematic uncertainties are incorporated into the TGC fit by allowing the expected num- 
bers of signal and background events to vary in the fit and constraining them to have a Gaussian 
distribution around their expected values with their systematic errors taken as the width of the 
distributions. The systematic uncertainties, excluding those on efficiency and background, are 
assumed to be correlated between the three different event selections. 



Data from lower centre-of-mass energies are included assuming all systematic errors to be fully 
correlated between energies. The corresponding log L curves are used in combination with the 
results of the angular distribution analyses which are described in the following sections. The 
full set of results are then presented in Figure |12|, where the cross section contributions are 
shown as dotted lines. 



5 Analysis of the W + W — > qqlv? angular distributions 

The analysis of the W + W~^ qq£i^ channel is performed in three different ways using optimal 
observables (OO), a binned maximum likelihood (BML) fit, and the spin density matrix (SDM). 
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These three methods will be described in the following sub-sections. All methods use the same 
selection and reconstruction procedure. The selection procedure starts from the event sample 
used for the total cross-section analysis with further cuts imposed in order to assure a reliable 
reconstruction of the event kinematics and to reduce further the background. 



5.1 Event selection and reconstruction 



The selection of qqeF e and qq/iX 7 ^ events for the total cross-section analysis results in a single 
track being identified as the most likely lepton candidate. The electron momentum vector 
is reconstructed by the tracking detectors and the energy is measured in the electromagnetic 
calorimeters. In the case of muons, the momentum measured using the tracking detectors is 
used. The remaining tracks and calorimeter clusters in the event are grouped into two jets 



using the Durham k± algorithm ||32|| . The total energy and momentum of each of the jets are 
calculated with the method described in [ |33|| . A one-constraint kinematic fit is then performed 
on the events, requiring energy-momentum conservation and allowing for a massless neutrino. 
Events are accepted if this fit converges with a probability larger than 0.001. This cut rejects 
about 2% of the signal events and 10% of the background. 

To improve the resolution in the angular observables used for the TGC analysis, we perform 
a further kinematic fit with three constraints requiring energy-momentum conservation and 
the fitted masses of both the hadronic and the leptonic systems to be constrained to the 
average Tevatron W mass, M w =80.40 GeV/c 2 |17| , within the W widthF]. We demand that 
the kinematic fit converges with a probability larger than 0.001. For the ~3% of events which 
fail at this point we revert to using the results of the fit without the W mass constraint. 

The selection of the qqrXv events results in the identification of the most likely tau decay 
candidate classified as an electron, muon, one-prong hadronic or three-prong hadronic decay. 
The remaining tracks and calorimeter clusters in the event are grouped into two jets as described 
above. However, these events cannot be reconstructed in the same way since there is more than 
one unobserved neutrino. Nevertheless, as the tau is highly relativistic, its flight direction can 
be approximated by the direction of its charged decay products, leaving us with four unknown 
quantities, the tau energy and the prompt neutrino three-momentum vector. A one-constraint 
kinematic fit is performed on the event requiring energy-momentum conservation and equality 
between the masses of the hadronic and leptonic systems. This fit is required to converge with 
a probability larger than 0.001. In addition, the angle between the two jets in the hadronic 
system, 9jj, is required to satisfy cos 9jj < —0.2, and the angle between the tau and the closest 
jet is required to exceed 20°. These cuts reject 20% of the signal and 42% of the background. 
Furthermore, these cuts suppress those qqrI7 r events which are correctly identified as belonging 
to this decay channel but where the tau decay products are not identified correctly, leading to 
an incorrect estimate of the tau flight direction or its charge. The fraction of such events in 
the qqrz/ T sample is reduced by these cuts from 13% to 8.5%. 

After these cuts, 324 qq?z^ candidates remain (135 qqez7 e , 116 qq/iX 7 ^ and 73 qqrz7 T ). The 

10 The W mass distribution is treated as a Gaussian in the kinematic fit. However, in order to simulate the 
expected Breit-Wigner form of the W mass spectrum, the variance of the Gaussian is updated at each iteration 
of the kinematic fit in such a way that the probabilities of observing the current fitted W mass are equal whether 
calculated using the Gaussian distribution or using a Breit-Wigner. 
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remaining background fraction, assuming Standard Model cross-sections for signal and back- 
ground processes, is 6.2%, not including cross-migration between the three lepton channels. 
The background sources are: four-fermion processes after subtracting the contribution of the 
CC03 diagrams (3.8%), Z°/ 7 -> qq (1.9%), W+W"-> qqqq and lv t l'V v (0.3%) and two-photon 
reactions (0.2%). 

In the reconstruction of the qq£z/^ events we obtain cos #w by summing the kinematically fitted 
four-momenta of the two jets. The decay angles of the leptonically decaying W are obtained 
from the charged lepton four-momentum, after boosting back to the parent W rest frame. In 
the hadronically decaying W we are left with a twofold ambiguity in assigning the jets to the 
quark and antiquark. This ambiguity is taken into account in the analyses described below. 

In Figure |3] (a-e) we show the distributions of all the five angles obtained from the combined 
qq£i/£ event sample, and the expected distributions for Ag\ = ±1 and 0. These expected 
distributions are obtained from fully simulated Monte Carlo event samples generated with 
Excalibur, normalised to the number of events observed in the data. Sensitivity to TGC 
is observed mainly for cos^w. The contribution of cos 6^, 0|, cos#* et and <j>* to the overall 
sensitivity enters mainly through their correlations with cos#w 



5.2 The Optimal Observable (OO) analysis 

The concept of optimal observables jlTj] is used to project the five kinematic variables of each 
event onto a single observable. It has been shown that for a differential cross-section which 
is linear in the parameter to be determined, the sensitivity of the optimal observable is the 
same as for a multi-dimensional maximum likelihood fit. The W-pair differential cross-section 
is, however, a second-order polynomial in the couplings, thus leading to a loss of sensitivity 
when constructing the optimal observable without the quadratic term. Nevertheless, in the case 
where the deviations of the couplings from zero are small, the quadratic terms are suppressed 
and the loss in sensitivity is far outweighed by the great advantage of a parametrisation in only 
one dimension. This method also allows usage of all the available kinematic information from 
the event, whereas in other methods (see the next sub-sections) some part of the information 
is lost due to computational limitations. 

The optimal observable for measuring any particular TGC parameter, a, is constructed for each 
event i with the set of phase-space variables fl = (cos6\y, cos0|, </>£, cos#* et , 0* t ) by differentiat- 
ing the differential cross-section for the event, crj(a) = da(fl, a)/dQ)\m = Q^, with respect to the 
parameter a evaluated at the Standard Model value, a—0, and normalising to the Standard 
Model cross-section, 

1 da,i(a) 

The differential cross-section used here is the Born cross-section, taking into account the twofold 
ambiguity in the definition of the measured cos 9* and (f>* et . 

Figure ^(f) shows the distribution of the optimal observable corresponding to Agf, along with 
the predictions corresponding to the Standard Model and to Ag\ = ±1, as obtained from corre- 
sponding Monte Carlo samples. Similar distributions are obtained for the optimal observables 
corresponding to Ak 7 and A. 
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A binned maximum likelihood fit of the expected optimal observable distributions to the data is 
performed to extract each coupling, assuming the other two couplings to have Standard Model 
values. The expected optimal observable distributions for various couplings are derived from 
Monte Carlo events and normalised to the number of events in the data in order to exclude any 
information from the overall production rate in this part of the analysis. To obtain the optimal 
observable distribution at intermediate parameter values, a reweighting technique is applied 
to sets of Monte Carlo events generated with different couplings using a Born-level differential 
cross-section but taking into account ISR effects. 

The method is tested with about 350 Standard Model and 200 non-Standard Model Monte Carlo 
subsamples, each corresponding to the collected luminosity. Summing up the log L curves from 
the different subsamples corresponding to the same coupling gives results which are consistent 
with the coupling value used in the Monte Carlo generation. Unfortunately, for the Standard 
Model subsamples the logL curves corresponding to Ak 7 tend to have two minima, one close 
to the generated value, Ak 7 =0, and the second around Ak 7 =2.5. The second minimum is 
due to the fact that, for example, the normalised differential cross-section for Ak 7 =2.5 is more 
similar to the Standard Model than the one for Ak 7 —1. In about 18% of the cases the wrong 
minimum is the deepest one. However, the cross-section information already excludes a value 
around 2.5 for Ak 7 (see Figure |1^). Therefore, the minimum at that point can be disregarded. 
A similar problem occurs in 3% (7%) of the Standard Model subsamples also for Agf (A) where 
the constraint from the cross-section information is even more stringent, so it is handled in the 
same way. 

The error interval is defined in the usual way, as the region where the logL function is higher 
than its minimum by no more than 0.5. To test the reliability of this error estimate, the 
fraction of subsamples where the correct value is within the error interval is calculated. The 
error estimate is considered to be reliable if the calculated fraction is consistent with 68%, 
which is the case for all couplings. 

Performing the fit to the 183 GeV data, we obtain the results quoted in the first row of Table || 
The statistical errors obtained from the fit are consistent with the expected values which are 
estimated from the distribution of the Monte Carlo subsample fit results. 

Biases introduced by uncertainties in the simulation of the detector acceptance are checked 
by removing events in phase-space regions where the efficiency changes rapidly. Cuts on the 
polar angles of the charged lepton and the neutrino, on the angle between the lepton and the 
nearest jet, as well as a cut on the lepton energy, were introduced. Furthermore the range of 
the OO-values that are considered in the fit is changed by approximately ±10%. The variations 
of the fit result due to these changes fall within the statistical fluctuations determined using 
many Monte Carlo subsamples with the size of the data sample. 

The effect of Monte Carlo statistics is studied by allowing the expected number of events in 
each bin of the optimal observable distribution to vary around its central value within its error 
according to a Gaussian distribution. This effect is found to be negligible. 

Additional systematic uncertainties are studied due to the following effects. 

a) For quark jets in the Monte Carlo samples, the resolutions of the three jet parameters 
(energy, cos^,^) are varied by 10%, and the energy scale is shifted by 0.5% to account 
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for systematic uncertainties in the modelling of the jet reconstruction. The size of the 
variations is determined from extensive studies of back-to-back jets at LEP Z° energies. A 
possible systematic shift in the reconstructed direction of the W boson has been estimated 
using radiative Z / 7 — > qq events. The shift in | cos #w | was found to be less than 0.01 Jll] . 
The same reconstruction uncertainties are assumed also for the r jets in qqrI7 r events. 

b) Uncertainties in the fragmentation model are studied using a Monte Carlo sample gener- 
ated with Herwig rather than Pythia. 

c) Differences between Monte Carlo generators not related to fragmentation are investigated 
by replacing the Excalibur reference sample by Koralw. This latter generator also 
uses the Jetset fragmentation scheme but has a different treatment of ISR and simulates 
correctly the r helicity effects. 

d) Systematic effects arising from the simulation of the Z°/7 — > qq background are taken 
into account by replacing Pythia with Herwig. The two-photon background is removed 
to test its impact on the measurement. The four-fermion processes from e + e~ — > Wez7 e , 
e + e~ — > Z°e + e~ and e + e~ — > Z°Z° are added as background to the CC03 diagrams, ne- 
glecting interference. The justification for this simplification is verified by using a fully 
simulated four-fermion Monte Carlo (Excalibur) as the reference sample. In addition, 
subsamples of fully simulated four-fermion Monte Carlo events have been fitted and no 
bias was found. To test the Monte Carlo generators, samples of Koralw and grc4f 
generated with the CC03 W-pair diagrams according to the Standard Model are used as 
test samples and all the TGC results are found to be consistent with zero. 





Source 


Error on parameter 






Ak 7 


Agl 


A 


a) 


Jet reconstruction 


.114 


.029 


.033 


b) 


Fragmentation 


.088 


.012 


.012 


c) 


MC generator 


.049 


.003 


.011 


d) 


Background 


.093 


.029 


.010 


Combined 


.178 


.043 


.038 



Table 7: Contributions to the systematic uncertainties in the determination of the TGC 
parameters in the OO analysis for the W + W~ — > qq£z^ angular distributions of the 183 GeV 
data. 

Table [7] summarises the various components of the systematic errors. These are incorporated 
into the logL functions by convolving their distributions, assumed to be Gaussian, with the 
likelihood functions. The corrected likelihood functions resulting from this convolution are used 
to obtain new values for the TGC parameters with modified errors which include the systematic 
uncertainties. The new TGC values may differ from those obtained before the convolution due 
to the asymmetric nature of the likelihood function. The results are listed in Table 



5.3 The Binned Maximum Likelihood (BML) analysis 

The BML analysis has been used to analyse our previous data taken at 161 GeV and 172 GeV 
and is fully described in || |TTJ. The finite statistics of the data and Monte Carlo samples 
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limits the number of kinematic observables which can be used to three. The most sensitive 
ones, namely cos#w> cos 9% and <$>\ are chosen^. The expected three-dimensional distribution 
of these observables as a function of the TGCs is fitted to the corresponding distribution of 
the data. The expected distribution is obtained from MC events generated according to the 
Standard Model with the CC03 diagrams. These events are reweighted to correspond to any 
other TGC value. The reweighting procedure uses very large statistics MC samples at the 
generator level, generated with ten different sets of TGC values and used to produce ten basis 
distributions. Since the differential cross-section is a second-order polynomial in the TGCs 
(see section [O] above), the distribution corresponding to any other set of TGC values can be 
obtained by a proper linear combination of the ten basis distributions. 

This procedure accounts for the effects of efficiency, resolution and background from other 
W + W~ decay channels. To account also for non-W + W~ background, the contribution of 
background sources which do not depend on the TGCs, namely Z°/7 — > qq, qqe + e~ and two- 
photon reactions, is added to the expected distribution used in the fit. The distributions from 
these sources are obtained from corresponding Monte Carlo samples. Background from other 
four-fermion diagrams which is TGC-dependent is neglected at this stage, and biases caused 
by this neglect are evaluated by performing the fit on four-fermion MC events taken as a test 
sample. These biases are subtracted from our fit results. 

The expected distribution is normalised to the number of events observed in the data in order 
not to incorporate any information from the overall production rate in this part of the analysis. 
The probability for observing the number of events seen in each bin is calculated using Poisson 
statistics. The statistical fluctuations in the Monte Carlo are taken into account using the 
method of reference | ^| . Consequently, the errors obtained from the fit include both data and 
Monte Carlo statistics. 

Detector level MC samples generated with different TGCs are used to verify that the BML 
fit method introduces no bias. The reliability of the statistical errors obtained from the fit is 
checked with many Monte Carlo subsamples, in the way described in section |5.2| . A problem 
occurs only for Ak 7 where the fraction of subsamples with the correct TGC value within the 
error interval falls below 68%. It is found that dividing the logL function by a correction factor 
of 1.21 increases the error interval in such a way that this fraction reaches 68%. This correction 
factor is then applied to the logL function extracted from the data. The statistical errors 
obtained in this way are also consistent with the expected values which are listed in Table |8|. 

The systematic errors are evaluated and incorporated into the log L functions in a similar way 
as for the 00 method. The results are consistent with those of the 00 method. 



5.4 The Spin Density Matrix (SDM) analysis 

Spin density matrix elements are observables directly related to the polarisation of the W 
bosons [0,0 • Additional insight into the underlying physics may be gained from these spin- 
related observables, and the relative production of various helicity states of the W bosons can 
be measured. Comparing the spin density matrix elements with the theoretical predictions 

11 The observables cos#* t and 0* t have smaller sensitivities due to their twofold ambiguity. 
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allows a model-independent test of the TGCs. If deviations were detected, this method would 
give information on the structure of the anomalous couplings. 



Spin density matrix elements are normalised products of the helicity amplitudes J 7 ^ (cos 9-w) , 
where r_ and r + are the helicity states of the W~ and W + boson, respectively, and A denotes the 
spin of the e + e~ system. The two-particle joint density matrix elements are defined according 
to 

Pr-r'_r + r' + (COS W ) = — ( A) • 

Z^At + T_ K T-T+ I 

Due to limited statistics in the analysis described here the single-W spin density matrix elements 
p TT i are used. The spin density matrix of the W~ boson, for example, is obtained from the 
two-particle joint density matrix elements by the relation 

P^_ T ,_(C0S9 W ) = J2Pt-t>_t + t+ (cos w ). 

The matrix p rr / is hermitian, thus having six independent matrix elements. The diagonal 
elements p TT of the spin density matrix are real and can be interpreted as the probability to 
produce a W boson with helicity r. The off-diagonal elements are complex in general, but for 
CP-conserving theories the imaginary parts vanish. 



The spin density matrix elements are extracted with the help of projection operators fl|.|43 



These operators reflect the standard V — A couplings of the fermions to the W boson in the 
W decay. The single-W spin density matrix elements p TT > can be extracted using the threefold 
differential cross-section da/dcos6^dcos6*d(j)* from the relation 



Wj/ Q ,da(e + e~ -> WiW 2 ) _ 1 f da(e + e~ -> W 2 ff) 



rfcos W £> r Wi -»ff J dcos fc'w d cos 0\ d(p* 

Here A Tr / is the suitable projection operator for extracting the spin density matrix element p TT i 
and Br Wl _»ff is the branching ratio for the decay of the W boson considered. Expressions for 
the projectors can be found in j43|. Experimentally, equation ^corresponds to 

1 N 

p Tr /(cos6» w ) = — 2jA TT ,(cos0*j,0*i) 
^ v i=i 

where N is the number of events and A = A(cos#*j, </>*J is the value of the projection operator 
of event i. The spin density matrix elements are extracted in bins of cos^w- 

For the leptonically decaying W, the full information about the W decay angles is accessible 
and the matrix elements can be directly extracted according to @. For hadronic decays of the 
W boson where the quark cannot be distinguished from the antiquark, only the folded angular 
distributions of cos 8* and <fi* are directly available. These folded distributions can be identified 
only with the symmetric^ part of the angular distributions M, whereas no information about 
the antisymmetric part is accessible. The following symmetric (combinations of) spin density 
matrix elements can be extracted from the folded angular distribution by using the symmetric 
part of the projection operators only: 

P++ + P-, Poo, Re(p+_), Im(p+_), Re(p +0 -p- ), Im(p +0 + P-o)- 



Symmetric/ antisymmetric under the transformation cos#* — > -cos( 
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The correlations between the single-W density matrix elements are calculated analytically. 

In order to compare the density matrix elements extracted from data with the theoretical 
predictions, the data have to be corrected for experimental effects, like selection efficiency, 
acceptance, angular resolution and the effect of ISR. A simple approach, based on a Standard 
Model Monte Carlo, is performed. First the expected background, taken from Monte Carlo, 
is subtracted on a statistical basis. Subsequently the data are corrected by multiplying the 
three-dimensional angular distribution by a correction function which is given by the ratio of 
the angular distributions for the reconstructed angles in the selected events and the generator 
level distributions of all events. The effect of ISR is approximately accounted for by reducing 
the centre-of-mass energy by the mean energy of initial state photons, (-Eisr) = 1.6 GeV, as 
determined from Monte Carlo. 

Distributions of the spin density matrix elements of the leptonically decaying W and the hadron- 
ically decaying W are shown in Figures || and |5|. Overlaid are the analytical predictions as ex- 
pected in the Standard Model. The errors are statistical. The systematic errors are estimated 
to be less than 10% of the statistical errors. This is verified by replacing the Standard Model 
Monte Carlo sample used for the correction of detector effects with samples generated with 
Agf = ±1 values which are far outside the allowed region obtained in our fit (see below). 

The measurements of the W-pair cross-section, the cos^w distribution and the spin den- 
sity matrices can be combined to give a measurement of the semi-inclusive differential cross- 
sections, to produce a transversely polarised W, e + e~— >WtW, or a longitudinally polarised 
W, e + e~— >WlW, where in either case the second W can have arbitrary helicity. For this pur- 
pose, the raw cos#w distribution, plotted in Figure ||(a), is corrected for detector effects by 
subtracting the background and multiplying by a correction function for efficiency and resolu- 
tion obtained from Standard Model Monte Carlo. The systematic error associated with this 
correction is estimated in the same way as for the 00 analysis. In addition, as was done 
for the spin density matrix elements, the Standard Model Monte Carlo sample used for the 
correction of detector effects is replaced by samples generated with Ag\ = ±1. The corrected 
cos #w distribution is multiplied by our measured total cross-section and the corresponding spin 
density matrix elements after combining their values from the leptonic and hadronic decays. 
The resulting differential cross-sections, plotted in Figure [], are seen to be consistent with the 
Standard Model expectations. Integrating over these cross-sections, the overall fraction of lon- 
gitudinally polarised W bosons is determined to be 0.242±0.091±0.023, where the systematic 
error is dominated by uncertainties in the jet and tau resolutions (0.017) and MC generator 
(0.015). The expected value for this fraction is 0.272 (0.392, 0.405) for the Standard Model 
(Atf=+1,-1). 

The spin density matrix elements can be used to extract the TGCs by comparing them, before 
the various corrections, to those expected from Monte Carlo events for different TGC values. 
This is done with a reweighting technique, taking into account all experimental effects. The W 
production angle and the spin density matrix elements of both W bosons are used in the fit. 
We obtain results consistent with those from the OO method. 
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5.5 Summary of the qqiug analyses 



All three methods described above give consistent results for the 183 GeV data. The 00 results 
are chosen to be combined with the results of the other channels and the cross-section result 
since this method is expected to give the best sensitivity, as can be seen from the expected 
errors listed in Table B 



Level of results 


Ak 7 




A 


Without systematics 


-n 91 +0 - 61 


0.00 ±0.14 


-0.19 ±0.15 


Expected errors, 00 method 
Expected errors, BML method 
Expected errors, SDM method 


±0.46 
±0.56 
±0.50 


±0.13 
±0.16 
±0.14 


±0.14 
±0.19 
±0.16 


Including systematics 


n 1 S+0-77 
-U.18_ a56 


0.00 ±0.16 


-0.19 ±0.16 


Including 161/172 GeV angular data 


-0.08185 


-0.01 ±0.15 


-0.16 ±0.16 


Including cross-section 




-0.01 ±0.14 


_n 14+ - 15 

u - iq: -0.14 



Table 8: Measured values of the TGC parameters using the qq_z/£ event selection analysed with 
the OO method, with and without taking into account the systematic uncertainty. We also list 
the results after combining with our 161 and 172 GeV data and after combining with the cross- 
section information of qqZz/^ data from all centre-of-mass energies. The expected statistical 
errors obtained from Monte Carlo studies are listed for all three analysis methods. 



Although our method of OO analysis using the optimal observable distribution is preferred for 
single parameter fits, it is less well suited to a fit of two or three TGC parameters. In this 
case, two- or three-dimensional distributions of the optimal observables corresponding to the 
different TGC combinations would be required, and it is simpler to fit the angular variables 
directly, using the BML method. 

The 183 GeV results are combined with our 161 GeV and 172 GeV results [O. In the 161 GeV 
analysis only the parameter aw<t> has been analysed. Therefore, this analysis is extended 
here using the same (BML) method to include the three parameters Ak 7 , Agl and A, along 
with the two-dimensional and three-dimensional fits. All the tests and the systematic studies 
are done in the same way as for the 183 GeV BML analysis. The correlation between the 
systematic errors of the qqZz/^ results from the three centre-of-mass energies is neglected since, 
using the BML method, it is found to affect the results by no more than 10% of their statistical 
errors. 

Table |8] summarises the qq£z/^ results from the 183 GeV data, the results obtained after combin- 
ing the 161 and 172 GeV data, and the combined qq.i^ results obtained after adding also the 
cross-section information for that channel from all centre-of-mass energies. The logL curves 
corresponding to the combined qq£f£ results are shown in Figure |TT[ The correlation between 
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the systematic errors of the angular distributions and the cross-section, mainly due to the un- 
certainty in the background level, is neglected as it affects the results by less than 1% of their 
statistical errors. 



6 Analysis of the W + W — > qqqq angular distributions 

The analysis of the W + W - -^ qqqq channel is performed on both the 172 GeV and the 183 GeV 



data. It uses the same event samples as the cross-section analyses (see section |3]3| and ref. JTTI ). 



6.1 Event selection and reconstruction 

Using the Durham k± algorithm | 3^| , each selected event is forced into 4 jets, whose energies 



are corrected for the double counting of charged track momenta and calorimeter energies [33 
To improve the resolution on the jet four-momentum we perform a kinematic fit requiring 
energy-momentum conservation and equality of the masses of the two W candidates (5-C fit). 
The event reconstruction is complicated by the ambiguity in the choice of the correct di-jet 
combination and uncertainties in the determination of the W charge. The latter is assigned by 
comparing the sum of the charges of the two jets coming from the same W candidate. Each W 
charge is defined as: 

n Ni + Nj N k + N t 

WW! — ^4 , WW 2 — ^4 



V 4 D V 4 D 

Z^m=l - Ly m Z^ m =l 

where jets i and j belong to Wi, jets k and I belong to W2, and Nj Gt and Dj et are given by: 

TV N 

N J et = J2 ( lM\\\ - 5 , D jet = J2\Pi\\\°- 5 

i=l i=l 

where is the charge of the i th track, p^u is the projection of its momentum along the jet axis 
and N is the total number of tracks in the jet. 

One di-jet combination is chosen out of the three possible ones by a likelihood algorithm. The 
input variables to the likelihood are the di-jet invariant masses, obtained by a kinematic fit 
requiring energy and momentum conservation (4-C fit), the charges of the two W candidates, 
and the probabilities of the 5-C kinematic fits. Once a jet combination is chosen, the W~ is 
defined as the di-jet whose charge is more negative. According to the Monte Carlo simulation, 
these choices correspond to a probability of selecting the correct di-jet combination of about 
78% (82%) at 183 (172) GeV, and to a probability of correct assignment of the W charge, 
once the correct pairing has been chosen, of about 76%. Both probabilities show a few percent 
dependence on the value of an eventual anomalous TGC. Incorrect jet pairing and wrong de- 
termination of the W charge decrease the sensitivity of the reconstructed cos^w distribution 
to possible TGCs. To increase the fraction of events correctly reconstructed, we apply fur- 
ther cuts on the value of the jet pairing likelihood, L > 0.8, and on the W charge separation, 
IQwi ~Qw 2 1 > 0.04. The distributions of these variables, for the data collected at 183 GeV and 
for the corresponding Monte Carlo, are shown in Figure [7]. The cut values are the result of a 
compromise between loss in efficiency and gain in sensitivity. After these cuts, the probability 
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of correct pairing increases to 86% (87%) at 183 (172) GeV and the probability of correct de- 
termination of the W charge to 82%. The selection efficiency is reduced by about 40%, whereas 
the Z°/7 — > qq background is reduced by 60% and 70% at 172 and 183 GeV respectively At 
172 GeV, 41 events survive all cuts out of the 59 initially selected, and 241 out of 438 events 
remain finally at 183 GeV. The cos 6^ distributions obtained at both energies are shown in 
Figure |8| together with the distributions predicted by the Monte Carlo for different values of A. 



6.2 TGC analysis 

Due to lack of separation between quarks and antiquarks, in W + W _ — > qqqq events most of the 
sensitivity to the TGCs is contained in the cos 6^ distribution. A binned maximum likelihood 
method is thus used in this channel. The measured cos 6*w distribution is divided into ten bins 
in the range [-1,1]- The expected cross-section in each bin i, di, due to W + W~ production is 
parametrised as a second-order polynomial in the analysed coupling. Monte Carlo reference 
histograms, including all effects of acceptance, resolution, incorrect jet pairing and incorrect 
determination of the W charge, are used to build the parametrisation. The contribution of 
the Z°/7 — > qq background as obtained from the Pythia Monte Carlo is added to the W-pair 
expectation. To avoid any dependence on the overall production rate, the expected cross-section 
is normalised to the total number of candidates. The expected distribution is fitted to the data 
using a binned maximum likelihood method and assuming Poisson statistics. The result of the 
fits to Ak 7 , Agl and A as obtained from the 183 GeV data sample are given in the first row of 



Table 10 where the quoted errors are statistical only. 



The reliability of the statistical error estimates is checked by performing the analysis on many 
subsamples of Monte Carlo events, each corresponding to the data luminosity. The r.m.s. 
values of the fit results corresponding to these subsamples are quoted in the second row of 



Table [10]. These values are somewhat lower than the statistical errors coming from the fits to 
the data. However, the data statistical errors are still compatible with their expected values 
with probabilities above 5%, as obtained from the MC subsample study. As a further test, the 
fraction of subsamples where the correct TGC value turns out to be inside the error interval is 
checked and found to be consistent with 68%, as required. 

The following sources of systematic errors are considered: 

a) The effect of possible differences in the jet resolution between data and Monte Carlo is 
studied in the Monte Carlo by smearing and shifting the jet energies and directions in the 
same way as done for the qqZi^ analysis. The resulting changes caused to the measured 
TGCs are added in quadrature and taken as a systematic error. 

b) Possible dependences on fragmentation models are studied by comparing the results of 
the fit to two Monte Carlo samples generated with Pythia. In the first sample, the 



Jetset [22] fragmentation model is implemented as for the reference Monte Carlo sam- 
ples, whereas in the second it is replaced by Herwig. The implementation of Herwig 
mainly results in a higher probability of correct jet pairing and correct W charge with 
respect to Jetset. The statistics of the data are insufficient to discriminate between the 
two models. We assign as a systematic error for each coupling the difference we find in 
the results of the fits to the two Monte Carlo samples. 
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c) We study possible biases due to the choice of the Monte Carlo generator by repeating the 
fit to W + W" samples generated with Pythia, Koralw and grc4f . 

d) A systematic error on the estimation of the background is determined by varying both its 
shape and normalisation. The shape predicted by Pythia is replaced by that predicted by 
Herwig. The normalisation is varied by the background uncertainty as determined from 
the analysis of the cross-section. The feed-through from other W boson decay channels 
is less than 0.2% and is neglected. The effect of neglecting the contribution of four- 
fermion diagrams is studied by performing the analysis on various four-fermion Monte 
Carlo samples generated with EXCALIBUR and grc4f . The corresponding shifts on the 
fitted values of each TGC parameter are added in quadrature and taken as a systematic 
error. 

e) Bose-Einstein correlations (BEC) between bosons originating from different W in the 
event might affect the measured W charge distribution. This effect is studied at 183 GeV 
with a Pythia Monte Carlo sample where BEC have been implemented. 



f) The measured W charge distribution might be affected also by colour reconnection. This 
effect is investigated at 183 GeV with Monte Carlo samples generated with Pythia, in 
two classes of reconnection models called type I and type II [fl5 |, and with Ariadne |[46[| . 



g) Finally, a possible presence of biases in our analysis or fit method is checked by analysing 
Monte Carlo events generated with known values of the couplings, ranging from -2 to 
+2. As a consistency check we use large Monte Carlo samples. We also perform fits to 
many small Monte Carlo subsamples as mentioned above and examine the sum of their 
likelihood functions. 





Source 


Error on parameter 










A 


a) 


Jet reconstruction 


0.26 


0.08 


0.09 


b) 


Fragmentation 


0.30 


0.10 


0.20 


c) 


MC generator 


0.45 


0.10 


0.13 


d) 


Background 


0.22 


0.06 


0.10 


e) 


BEC 


0.07 


0.04 


0.02 


f) 


Colour reconnection 


0.01 


0.01 


0.01 


g) 


Fit bias tests 


0.20 


0.08 


0.09 


Combined 


0.68 


0.19 


0.29 



Table 9: Contributions to the systematic errors in the determination of the three TGC 
parameters for the W + W~— > qqqq angular distribution analysis of the 183 GeV data. 



The results of all these systematic studies for the 183 GeV data are summarised in Table |9|. The 
systematic errors for the 172 GeV data are similar and are assumed to be fully correlated with 
those of the 183 GeV data. The combined systematic errors are convolved with the likelihood 
functions. The results for the 183 GeV data sample after convolving the systematic errors 
are listed in the third row of Table 10|. After combining with the 172 GeV data, we obtain 
the results listed in the fourth row of Table [10]. Finally, we combine also the results of the 
cross-section analysis in this channel, using the 161, 172 and 183 GeV data samples, yielding 
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the results quoted in the last row of Table [H]. The corresponding log L functions are plotted 
in Figure [TTJ. 



Level of results 


Ak 7 




A 


Without systematics 
Expected errors 
Including systematics 
Including 172 GeV data 
Including cross-section 


±1.13 

- l -- lD -1.34 

i-i5ll:| 

0.851SSS 


0-681SS 

±0.43 
n 70+0.88 

u - ' <J-0.72 

fi2+°- 97 

u - DZ -0.62 

24+ - 24 


7fi+ - 78 
u - ' °-0.67 

±0.52 
0.79l°;i 
0-681^ 

o.3ol°:l 



Table 10: Measured values of the three TGC parameters using the qqqq event selection from the 
183 GeV data, with and without taking into account the systematic uncertainties. The second 
row gives the expected statistical error from Monte Carlo studies. The fourth row shows the 
results after combining with the 172 GeV data, and in the last row we list the results after 
combining with the cross-section information of qqqq data from the 161, 172 and 183 GeV 
centre-of-mass energies. 



7 Analysis of the W + W — > ivtt'vp angular distributions 

In a W + W~^ Ivil'Vti event, there are at least two undetected neutrinos. Fortunately, in the 
small-width, no-ISR approximation, the W production and decay angles can still be recon- 
structed from the measured charged lepton momenta, provided there are only two missing 
neutrinos - thus restricting t =e or /i for both leptons. 

7.1 Event selection and reconstruction 

The selection procedure starts from the W + W~— > Jugi'Ve' event sample used in the cross-section 
analysis. We then apply further cuts, mainly to suppress the contribution of events with one 
or two r leptons which cannot be reconstructed. Only the 183 GeV sample is analysed, as the 
statistics in the lower energy samples are too small. 

Simple cuts are applied to ensure that the momenta and charges of the two leptons can be 
determined. We require the events to have at least two reconstructed cones containing charged 
particle tracks passing the quality cuts used in the primary Iv^'vp selection. The two highest- 
energy cones with charged particle tracks are assigned to the leptons and the charge of each 
cone is found by summing the charges of the tracks in the cone. The charges of the cones are 
required to be of different signs, or if one cone has a zero charge, the other must have a non-zero 
charge. 
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At this stage, 52% of the selected sample consists of events with at least one W— >rz/ decay. 
This fraction is suppressed by the following selections: 



The highest track momentum, pi, in each lepton cone must exceed 23.0 GeV, which is 
less by approximately twice the experimental resolution than the minimum momentum 
allowed for a lepton from an on-shell W decay in an event with no ISR at ^=183 GeV. 



The two lepton cones must each have no more than two good electromagnetic clusters [47 



Each of the lepton cones must be classified as an electron or a muon candidate, using the 
momentum of the most energetic track in the cone (pi), the electromagnetic calorimeter 
energy in the cone (Eem), and the hadron calorimeter energy in the cone (Eh c)> as 
follows: 

— electron candidate cones are required to have E EM /pi > 0.75 and E HC /pi < 0.1, 

— muon candidate cones are required to have E EM /pi < 0.1 and E HC /p\ < 0.5. 



After applying these r rejection cuts, the contamination from events with W— >ti> decays falls 
to 10%. This fraction is found to be essentially independent of the TGC parameters. 

In the data, 78 events are selected by the primary Iv^'Vy selection. Of these, 74 events pass 
the simple charge /momentum reconstruction requirements, and 30 pass also the r rejection 
selections. The purity of events with two charged leptons (e, fi or r) and two neutrinos in 
this final sample is predicted to be above 99%, where the main background (0.13±0.07 events) 
comes from the reaction e + e~—>e + e~ii + iJ,~. This level of background is negligible and is not 
considered further. 

The lepton momenta are reconstructed using the track momentum for muon candidate cones, 
and the energy seen in the electromagnetic calorimeter for electron candidate cones. 

The five characteristic angles are reconstructed event by event, in the approximation of zero 
W width and no ISR. This means that the two Eu systems are taken to have the W mass and 
the two W systems are assumed to recoil back-to-back in the laboratory frame, with a total 
energy equal to the centre-of-mass energy. When these constraints are applied, the six observed 
momentum components of the two charged leptons can be transformed into the five angles (plus 
an overall azimuthal angle which is of no interest). However, this transformation requires the 
solution of a quadratic equation, and thus there are either two, or no, real solutions for any 
specific event. The angle set {cos 6^, <fi* and (j)^} suffers from this ambiguity, but 9* and Q\ 
can be determined from the magnitudes of the lepton momenta alone. In the ideal case, where 
the W bosons are produced on-shell, and where there are no ISR or detector resolution effects, 
all events have two solutions: the ambiguity corresponds to a reflection ambiguity for the two 
neutrinos in the plane defined by the two charged lepton momentum vectors. The effects of W 
width, ISR and detector smearing can move leptons to momentum configurations not allowed in 
on-shell W decays yielding complex solutions for the momenta when attempting to reconstruct 
the neutrinos. These events lie preferentially in certain regions of the five-angle space, and 
losing them from the analysis would introduce biases. They can be recovered simply by taking 
the nearest real solution in five-angle space, taken to be the complex neutrino momentum 
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solution with the imaginary part set to zero. In these cases there is exactly one nearest real 
solution: the two solutions with complex momenta are complex conjugates of each other. 

Of the 30 selected candidates, 21 have two solutions found and the other 9 have their angular 
information recovered using the nearest real solution prescription. This is consistent with the 
Monte Carlo expectation of 8.9 events failing to have two reconstructed solutions. 



7.2 TGC analysis 



The fit for a TGC parameter a consists of minimising logL, defined as: 

/ dcr i 



log L = - lo § 



dn 

dcr 



a 



/£(a)/(fi)dn, 



where Q represents the five-angle set (equivalent to phase space for on-shell W bosons and no 
ISR) and / represents the acceptance, assumed to be independent of a (valid in the on-shell, 
no-ISR case). The differential cross-section da/dQ is determined using the program of Bilenky 
et al. ||. For events with two solutions for the five angles, the average of the differential 
cross-sections at the two solutions is used. 

The advantage of this unbinned maximum likelihood approach over binned techniques, as em- 
ployed in the qqqq and qqZz/g channel analyses, is that it uses the full five-angle information 
without any loss of information from binning. The main disadvantage is that the naive cross- 
section calculation of do j dVL does not include many of the effects which really occur: specifically 
the effects of ISR and W width, detector resolution, and backgrounds. The detector and selec- 
tion acceptance is partially included through the acceptance function f(fl) via an approximate 
analytical acceptance model. 

It is essential to demonstrate that this simple fit is effective in extracting couplings and that 
the effects not modelled in the fit function have a relatively small influence. The implications 
of unmodelled effects can be pernicious: they can give rise to biases in the fitted couplings and 
in the estimations of the fit errors. These issues are addressed via Monte Carlo tests, with 
high statistics to measure biases in the fit, and with large numbers of simulated low-statistics 
experiments to calibrate the errors. 

Studies are made with high-statistics Monte Carlo samples to test the biases in the fit method. 
These biases can arise from a variety of causes and are examined step-by-step by considering 
Standard Model samples with different effects included. Biases and systematic errors are derived 



as described in the following, with results shown in Table [IT]. In all cases the systematic errors 
are chosen so as to cover the differences in Monte Carlo models and the full size of any bias 
observed. In most tests, events from different Monte Carlo generators gave consistent results 
and the exceptions are noted. 



a) The effect of the primary Ipgi'Vy selection cuts is assessed by fitting with the true momenta 
of decay leptons, both before and after the selection cuts are applied, to fully simulated 
events. This test is made with Monte Carlo W-pair events generated by Excalibur, 
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grc4f and Koralw. The modelling of the main acceptance effects in the function / 
reduces the bias from the acceptance. 



b) Fits are made using the same fully simulated events, but fitting both with the true and 
the reconstructed momenta. 

c) The effect of the additional W— >rz/ rejection cuts on the Ttve&'Vgi signal (£ = e or /i) is 
assessed by fitting to the reconstructed simulated signal events before and after the cuts 
are applied. 

d) The W~— >rI7 T background is next included, and the simulated samples refitted. The size 
of the bias, averaged over the Excalibur, grc4f and Koralw samples, was found to 
be small. In this case, however, the different generators do give differences which are 
on the edge of significance: Excalibur and grc4f indicate a small positive bias to the 
results, Koralw a slightly larger negative one. The systematic error shown in Table [11] 
was chosen to cover all the different predicted biases. The modelling of r polarisation 
and decays is one area in which some of the generators are known to be defective, so that 
the different biases seen are not surprising: nonetheless we quote an error large enough 
to cover all three models. 

e) Four-fermion effects were tested by fitting Excalibur and grc4f samples produced either 
including only the CC03 W-pair diagrams, or when all four-fermion diagrams contributing 
to Iv^Vy final states were employed, including the contributions from extra final states 
such as z/ e z/ e . While the two Monte Carlo generators predict quite different changes 
in accepted cross-sections when the four-fermion effects are included, the shifts in fit 
results are small and compatible with each other. The central value of the bias was taken 



to be that of Excalibur, and is listed in Table O. The larger of the magnitude of the 
bias, and the difference between the shifts seen with Excalibur and grc4f , is taken as 
the associated systematic error. 



Source 


Ak 7 




A 


a) 


tvit'vii selection 


-0.03±0.12 


+0.05±0.05 


+0.01±0.05 


b) 


Detector resolution 


-0.07±0.07 


-0.05±0.05 


-0.02±0.02 


c) 


r rejection 


-0.07±0.07 


-0.05±0.05 


-0.02±0.02 


d) 


r background 


+0.02±0.04 


-0.01±0.07 


-0.02±0.03 


e) 


4-fermion effects 


-0.10±0.10 


-0.01±0.02 


-0.01±0.02 


f) 


Low statistics tests 


±0.56 


±0.05 


±0.03 


Combined 


-0.25±0.61 


-0.07±0.14 


-0.06±0.09 


logL correction factor 


1.17 


2.07 


1.39 



Table 11: Summary of biases, systematic errors and likelihood scale factors derived from 
Monte Carlo tests for fits to the lv$Ui> angular distribution. Rows a) to f ) show biases derived 
from the different tests discussed in the text, together with the associated systematic errors. 
The last row shows the log-likelihood scale correction applied. 

In addition to these detailed studies with Standard Model events, high-statistics bias tests were 
also made for non-standard couplings, to ensure that the fit biases are not changing strongly 
with the underlying TGCs. These tests were done separately for W-pair and full four-fermion 
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samples. In most cases, the best fits were obtained with values consistent with the generated 
couplings. Some fits, particularly for Ak 7 , show an additional minimum separated from the 
correct value, and in some cases this may fit slightly better than the generated coupling value. 
In one case (A/c 7 =l with the full four-fermion diagrams) the two minima merge. In all cases 
the true generated coupling was excluded by an amount equivalent to a log L interval of less 
than 0.25, with the present statistics. Incorrect second minima would in any case be removed in 
the full analysis by the addition of the event rate information and that from the other channels. 

The high-statistics studies described above evaluate the biases intrinsic to the unbinned max- 
imum likelihood fit method. It is also essential for this fitting approach to calibrate the fit 
errors, as they cannot be expected to be estimated free of bias. This is done in the same way 
as for the other decay channels, using many Monte Carlo subsamples with the same statistics 
as the data (30 event). A scale correction to the logL function is calculated such that in 68% 
of the subsamples the correct TGC value is inside the error interval. The scale corrections are 
listed in Table O 



The mean of the central values of these low-statistics fits were also compared, for each model, 
with the expected bias from the high-statistics test. The results are consistent for Ac/f and A 
but differ by 0.56 for Ak 7 . The larger of the differences or the statistical precision of the test 



is quoted as the "low statistics tests" systematic error, f), in Table [TT 



7.3 Results for the Iv^vy channel 



The reconstructed cos #w distribution is shown in Figure |9| and compared with the Standard 
Model expectation and different TGC hypotheses. The statistics is low, but the shape is con- 
sistent with the Standard Model expectation, as is the joint distribution of cos^ and cos Q\ 
shown in Figure 10. In both figures there is a clear difference visible between the different 
high-statistics Monte Carlo predictions. It is interesting to note that the correlation between 
the two 9* angles changes with A. Such an effect can only be measured well in the Iv^l'Vp 
channel because of the difficulty in distinguishing the quark from the antiquark direction when 
a W boson decays hadronically. 



Level of results 


Ak 7 




A 


Without systematics 
Expected errors 
Including systematics 
Including cross-section 


-1 00+ ' 94 

i - uu -0.92 

±1.08 
-0.74+1;?! 
-0.30^?° 


-0.791°;^ 

±0.71 
-0-761°;™ 
-0.29i°;ff 


-0.29+°^ 
±0.41 

n OQ+0.35 
U - ZO -0.39 

fl 1 fi+0-32 



Table 12: Fit results at different levels, in the iv^i'v^i channel. 



The results of the fits to the 30 selected candidates are listed in the first row of Table |T^. 
These results do not include any systematic errors but they correspond to the logL curves 
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after adjustment with the scale factor to obtain the correct statistical errors. The expected 
errors, defined to be the r.m.s. values of the Monte Carlo subsample results, are listed in the 



second row of Table [L2| and they are comparable with the statistical errors obtained from the 



corrected log L functions. 

In the next step, the logL functions are shifted to correct for the biases and convolved with the 
combined systematic uncertainties. The results are listed in the third row of Table [12|. Finally, 
the results are combined with those of the cross-section analysis corresponding to the tviH!Vy 
channel and using the 161, 172 and 183 GeV data samples. This combination is done by adding 



the corresponding log-likelihood curves, and the results are listed in the last row of Table 12. 



8 Combined TGC results 



The TGC results for the three event selections using the angular distributions are combined 
by summing the corresponding log L functions. The correlation between the systematic errors 
of the three results is neglected, since most of the important sources of systematic errors are 
relevant to a particular result and not common to all three of them. The logL functions 
obtained for the different couplings are shown in Figure [1^. Adding these functions to those 
obtained using the total cross-section yields the combined logL functions which are plotted in 
Figure [L^. Table [13] lists the combined results. 





Ak 7 


Ag\ 


A 


Combined results 
95% C.L. limits 


n 1 1+ - 52 

[-0.55, 1.28] 


o m+ 013 

U - Ui -0.12 

[-0.23, 0.26] 


-o.io±g : g 

[-0.33, 0.16] 



Table 13: Combined results of the three TGC parameters. The result on each parameter is 
obtained setting the other two parameters to zero. 



To study correlations between the three TGC parameters we also extract the logL as a function 
of all three variables, Ak 7 , Ag\ and A. Figure [13] shows the 95% C.L. contour plots obtained 



from two-dimensional fits, where the third parameter is fixed at its Standard Model value of 
zero. We also perform a three-dimensional fit, where all three couplings are allowed to vary 



simultaneously. The corresponding projections are plotted as dashed contour lines in Figure [13. 
As can be seen, the allowed range for each parameter is extended when the constraints on the 
other two parameters are removed. 

These three-dimensional fits can be used to obtain results for any other set of TGC parameters, 
as long as relations ([!]) between the five couplings are satisfied. For example, we determine the 
a parameters used in our previous publications to be a^=0. llloii?? a w<p = ~ 0-04 ± 0.09 

and aw=— OTOlo^f. The result on each a parameter is obtained assuming that the other two 
parameters vanish. 
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9 Summary and conclusions 



Using a sample of 877 W + W candidates collected at LEP at a centre-of-mass energy of 
183 GeV, we measure the total cross-section of, 

a(e+e-^W + W~) = 15.43 ± 0.61 ± 0.26 pb 

under the assumption that the W boson decay branching fractions and distributions of produc- 
tion and decay angles are all according to the Standard Model expectations. 

A measurement of the W branching fractions is made using the combined 161 GeV, 172 GeV 
and 183 GeV data samples. The results for the different leptonic decay channels are consistent 
with each other, as expected from lepton universality. Assuming Standard Model total cross- 
section and lepton universality, we obtain the hadronic decay fraction to be (67.9 ±1.2 ±0.5)%. 
From this result, a value for the CKM mixing matrix elements |V CS | = 0.99 ± 0.06 ± 0.02 is 
extracted, using also the measurements of the other matrix elements not involving the top 
quark. 

The total cross-section measurement, being consistent with the Standard Model prediction of 
15.72 pb -1 , can be used to place limits on anomalous triple gauge boson couplings. Those 
couplings are also investigated in an independent way, using the W-pair production and decay 
angular distributions for qq£z/^, qqqq and Iv^l'Vp final states. The W + W~^ qq?f£, being the 
most sensitive channel for this study, is analysed by three different methods, one of them 
utilising for the first time the spin density matrix of the W decay. All three methods lead to 
consistent results. The spin density matrix is used to measure the cross-sections to produce 
transversely and longitudinally polarised W bosons. Integrating over all angles, the fraction of 
longitudinally polarised W bosons is determined to be 0.242±0.091±0.023. 

The TGC measurements for all decay channels are combined, and the results obtained are, 

Ak 7 = ±0.111°$, 
Agl = ±0.0li°;{3, 
A = -0.10l°;& 

where each parameter is determined setting the other two parameters to zero. The precision of 
these results is comparable to the latest DO results obtained from boson pair production at the 
Tevatron pp collider ||. These results supersede those from our previous publications [p|,[TT 
They are all consistent with the Standard Model value of zero. 
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Figure 1: The distribution of the likelihood discriminant used to select events in the 
W + W~-^ qqqq selection is shown for all preselected events. The points indicate the data 
and the histogram represents the Monte Carlo expectation where the hatched area shows the 
estimated contribution of the total background. The selection cut is indicated by the arrow. 
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Figure 2: The dependence of <7ww on y/s, as predicted by Gentle forMyj = 80.40 GeV/c 2 . The 
W + W~ cross-sections measured at \fs = 182.7 GeV (this work), at \/s = 161.3 GeV jj^j and 
^/i = 172.1 GeV j2£j/ are shown. The error bars include statistical and systematic contributions. 
The dashed curve shows the expected cross-section if the ZWW couplings are zero. 
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Figure 3: Distributions of the kinematic variables cos #w, cos 6}, cos 9* ct , (/>%, 0* t and the optimal 
observable corresponding to Agf, as obtained from the qq£z/^ events. The solid points represent 
the data. The histograms show the expectation of the Standard Model (solid line) and the 
cases of Ag\=+1 and -1 (dotted and dashed lines respectively). The shaded histogram shows 
the non-qc(Iise background. Notes: 1. In the case ofW + ^iv decays the value of 4>} is shifted by 
7r in order to overlay W + and W~ distributions in the same plot. 2. The jet with < 0* ot < 7r 
is arbitrarily chosen as the quark (antiquark) jet from the decay of the W~ (W + j. 
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Figure 4: Spin density matrix elements for the leptonically decaying W as a function of cos #w 
The data points are corrected for experimental effects. The solid (dotted, dashed) lines show 
the predictions of models with Ag^=0 (+1, -I). All other anomalous couplings are set to zero. 
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Figure 5: Spin density matrix elements for the hadronically decaying W as a function of cos #w 
The data points are corrected for experimental effects. The solid (dotted, dashed) lines show 
the predictions of models with Agl=0 (+1, -I). All other anomalous couplings are set to zero. 
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Figure 6: Differential cross-section to produce a) a transversely polarised W and b) a longi- 
tudinally polarised W in a W-pair event where the second W can have any polarisation. The 
points represent the data and the solid (dotted, dashed) lines show the predictions of models 
with Agl=0 (+1, -1). The error bars include statistical and systematic uncertainties, except 
for a normalisation error of 4.3% associated with the total cross-section measurement. 
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Figure 7: For the qqqq channel, a) distribution of the output of the jet-pairing likelihood 
corresponding to the most likely combination for the data collected at 183 GeV (points) and for 
the Monte Carlo (histogram). The hatched area shows the contribution of wrong pairing, and 
the dark area represents the contribution of the background. The arrow indicates the cut value, 
b) Distribution of the charge separation of the W candidates for the data collected at 183 GeV 
(points) and for the Monte Carlo (histogram). The hatched area shows the contribution of 
correct pairing and incorrect W charge, the double hatched area shows the contribution of 
wrong pairing, and the dark area represents the contribution of the background. The arrow 
indicates the cut value. 
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Figure 8: Distribution of cos #w in qqqq events as obtained at 172 GeV (left) and 183 GeV 
(right). The solid points are the data. The Monte Carlo predictions for X=0,+2,-2 and 
X—0,+1,-1 are shown at 172 and 183 GeV, respectively, as solid, dashed and dotted lines. 
The hatched histograms show the contributions of the Z°/7 — > qq background, as predicted by 
Pythia. 
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Figure 9: Observed cos6> w distribution in the Iv^i'v^ analysis. All events enter with a total 
weight of one, but in the case of events with two ambiguous solutions for cos6>w, each solution 
enters with a weight of 0.5. 
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Figure 10: Observed distribution of cos 0* vs. cos 9% in the £vi£'V£> analysis, compared to 
the expected distributions from samples with standard (SM) and non-standard couplings, as 
indicated. In the panels corresponding to the Monte Carlo expectations, the area of each box 
is proportional to the fraction of the events falling into that bin. 
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Figure 11: Negative log-likelihood curves obtained in the different event selections: qq£z^ 
(dashed lines), qqqq (dotted lines) and Iv^'Vy (dash-dotted lines). Each curve is obtained 
by combining the results from the angular distributions and the cross-section and setting the 
other two TGC parameters to zero. Systematic errors are included. The solid line is obtained 
by combining the three event selections. 
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Figure 12: Negative log-likelihood curves obtained using different sources of information on the 
TGCs. The curves for each TGC parameter are obtained setting the other two parameters to 
zero. The dashed lines are obtained from the angular distributions, and the dotted lines from 
the total cross-section. All W decay channels are used and systematic errors are included. The 
solid line is obtained by combining the two sources of information. 
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Figure 13: The 95% C.L. two-dimensional correlation contours for different pairs of TGC pa- 
rameters. The solid lines are obtained by varying two parameters and fixing the third one to 
zero, which is the Standard Model value. The dashed lines show the projections of the three- 
dimensional conhdence regions obtained by varying all three parameters. The star indicates 
the best three-parameter fit values. These results are obtained from all cross-section data as 
well as angular distributions of all qqfz^ and 183 GeV qqqq data. 
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